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ABSTRACT 

A  formalism  is  developed  in  which  the  nuclear  potential  is 
explicitly  worked  into  the  radiation-matrix-elements.  The  result 
is  an  expansion  of  the  ordinary  matrix  element  in  a  power  series 
of  ^  and  In  each  term  of  this  series  occurs  explicitly 

pc 

the  potential  between  the  scattered  particles.  With  these  ex- 

4 

pressions  the  bremsstrahlung  of  n-He  -scattering  is  calculated 
if  the  nuclear  potential  is  assumed  to  be  a  square  well.  Further¬ 
more  the  bremsstrahlung  of  n-p-scattering  is  studied  for  two  kinds 
of  interaction;  a)  square  well,  b)  square  well  plus  hard  core. 

If  the  parameters  of  these  potentials  are  chosen  so  that  they  give 
the  same  scattering  length  and  effective  range,  there  remains  only 
a  small  difference  in  the  bremsstrahlung  cross  section  for  the  two 
types  of  potentials  (some  few  percent).  The  shapes  of  the  brems¬ 
strahlung  spectrum  are  equal  in  both  cases.  Some  structure  in  the 

shape  occurs  at  the  virgual  bound  state  of  the  s-wave.  The  magni- 

-30  -29  2 

tude  of  the  bremsstrahlung  cross  section  is  the  order  10  -  10  cm 

if  the  energy  resolution  is  0.2  MeV.  So  it  seems  that  nuclear  brems¬ 
strahlung  can  be  measured. 

♦Permanent  Address;  Physikalisches  Institut  der  University 

Freiburg/Brsg. ,  Freiburg/Germany. 


INTRODUCTION 


There  have  been  several  investigations  on  nuclear  bremsstrahlung 
1  -3 ) 

in  recent  years  .  The.  main  point  of  view  of  these  authors  was 
low  energy  bremsstrahlung  and  the  possibility  of  measuring  ^delay- 
times"  of  the  scattered  particles  due  to  compound  reactions.  Our 
point  of  view  on  nuclear  bremsstrahlung  is  different:  We.  ask 
whether  the  bremsstrahlung  cross  section  is  sensitive  to  certain 
features  of  the  nuclear  forces.  In  this  paper  we  especially 
investigate  the  differences  in  the  nuclear  bresstrahlung  cross 
section  due  to  a  hard-core-potential  on  the  one  hand  and  a  square- 
well-potential  on  the  other  hand.  These  potentials  are  taken  in 
such  a  way,  that  they  fit  the  same  low-energy  scattering  data; 
the  same  scattering  length  and  effective  range.  The  radius  of 
the  hard-core  is  taken  from  high-energy  scattering  experiments. 

In  order  to  investigate  this  question  we  develop  in  section  I 

a  method  so  that  the  nuclear  potential  is  explicity  worked  into 

the  radiation-matrix  elements.  What  we  get  for  the  usual  radiation- 

matrix  element  is  an  expansion  in  a  power  series  of  the  parameters 

—  and  ^2  where  p  is  the  nuclear  mass  and  hoo  the  photon¬ 
ic 

energy.  This  series  converges  very  rapidly  in  all  low-  and 
medium-energy  processes.  The  point  is  that  each  term  of  this 
series  contains  explicitly  the  nuclear  potential.  Therefore  a 
multipole  expansion  in  the  usual  way  is  possible  and  meaningful 
because  the  nuclear  potential  restricts  the  contribution  of  the 
matrix  element  to  nuclear  dimensions.  This  is  what  one  expects, 
that  even  in  scattering  processes  the  radiation  occurs  only  near 

-  2  -  • 
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the  scattering  center.  This  mathematical  method  i6  completely 

4  5) 

equivalent  to  the  well-known  Bloch-Nordsieck  ’  '  -transformation, 
but  we  put  it  in  a  lucid  light,  so  that  it  seems  that  insight 
is  deepened.  In  the  second  part  of  the  paper  we  discuss  the 
bremsstrahlung  of  a  square-well-  and  a  square-well-plus  hard-core 
potential  for  neutron-proton  scattering.  The  numerical  results 
are  done  only  for  s-scattering ,  but  it  is  straightforward  and 
only  a  numerical  task  to  calculate  the  contributions  of  other 
partial  waves.  The  expected  result  is  that  the  cross  section 
seems  to  be  different  by  few  percent  in  both  cases  (the  parame¬ 
ters  are  taken  so  that  the  scattering  length  and  effective  range 

are  equal).  On  the  other  hand  the  cross  section  is  of  the  order 
-3  -2 

10  -  10  millibarn  if  the  energy  resolution  is  taken  to  be 

0.2  MeV  and  therefore  measurable.  So  the  nuclear  bremsstrahlung 
seems  to  give  no  new  information  about  nuclear  potentials.  However 
other  potentials  which  are  spin  and  velocity  dependent  should  be 
discussed  and  their  effect  on  nuclear  bremsstrahlung  should 
be  studied. 


I. 

II:  The  Interaction  with  the  Radiation  Field. 


l/?e  .  firs  tn  consider. '.'the.. general  tint  e  taction  <&£etheotwo?:particles 
with  the  radiation  field  and  develop  a  form  of  the  interaction 
operator  which  is  suitable  for  practical  calulations.  We  have 


Hint 


* 


+ 


+  'y+jZ  (  ?x£)  c 


(l) 
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where  "p^  and  ^  are  the|  momenta;  >  mg  the  masses;  »  pg 

the  magnetic  moments;  ^  ,  7g  the  spin  operators  and  ,  ^g 

the  positions  of  the  two  particles.  Z  is  the  number  of  protons 

and  A  the  number  of  nucleons  in  the  recoil-nucleus.  T?  is  the 

wave  vector  and  T  the  polarization  vector  of  the  photon.  The 

first  term  in  (1)  is  the  *  /t  (£  is  the  vector  potential) 

interaction  of  the  two  particles  while  the  second  and  third  terms 

represent  the  Pauli-interaction  with  the  radiation  field.  The 

transformation*^  (1)  to  relative  &  "^2^  anc*  center-of- 

-  mii/ 1  ^  nuro 

mass  coordinates  (K  =  — — — r— -  )  gives 


Hint  -  -  |  ^ 


$  lwe 


7 rr/  _  niz  e 
M 


e 


- 1  +  2  e  ^  J  7T  e  * 


(2) 


+  iA*,o?-(5 


=*  I* 


*K)e 


it-(jph  -iff 


+F**?i(t‘k)e  J 


m^ffig 


Where  M  =  a,  ♦  m«  ,  it  =  - — -r—-  and  p,  P  are  the  momenta  of 
i  2  *  m^  ▼  mg 

the  relative  and  center-of-mass  motion  respectively.  The  matrix 
elements  of  the  second  term  in  (2)  are  small  of  the  order 
compared  to  the  first  one*  as  can  easily  be  seen  by  writing 
P/M  r  k  ;„a  _  where  "a*  is  a  length  of  nuclear  dimensions, 

*  8  C 

vCM  the  velocity  of  the  center-of-mass  and  hw  the  photonen:egrgy. 
This  term  represents  an  interaction  of  the  center-of-mass  motion 
with  the  ■internal*  relative  motion  via  the  radiation  field. 
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The  center-of-mass  coefficients  e 


-iff 


give  by  integration 


over  R-space  just  the  conservation  of  momentum,  except  for  the- 


second  term  in  {$)  where  we  get  the  additional  factor 


TT 


hj^wif)  is  the  momentum  of  the  center-of-mass  finally.  Usually 
when  calculating  nuclear  radiation  transition  probabilities  this 
term  can  be  neglected,  since  the  center  of  mass  is  nearly  at  tOst. 

We  shall  in  further  developments  drop  these  center-of-mass  dependent 


terms  and  use  the  interaction  in?-space: 

Hint  -  -  |  M  e 

c  u,  l  M 


.  e  tblti ■?  +2ear*r't'J]  + 

C  ft  L 


(3) 


+  i[r,*7  ■(?*?)£ 


-ck-Cjf 


?) 


+  rj.5?.(?*ic)e 
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In  calculating  transition  probabilities  for  a  nucleus,  the  bound 
state  wave  functions  differ  from  zero  only  in  a  region  raH  of 
nuclear  dimensions  and  one  has  k  •  a  «  1  .  Therefore  in  this 
Case  a  multipole  expansion  is  meaningful;  i.e.,  only  the  dipole 
and  quadriipole  transitions  occur.  That  this  is.  still  true  in 

i  ’ 

the  case  of  bremsstrahlung,  where  both  initial  and  final  wave 
functions  are  extended  over  the  whole  space,'  is  physically 
expected  but  not  easy  to  see  in  the  formula.  We  show  this  in 
developing  an  expression  for  the  matrix  element  C^lHint . 

4\ 

which  is  essentially  equivalent  to  the  well  known  Bloch-Nordsieck  )- 
transformation,  but  more  systematic  and  lucid. 
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Let  us  consider  for  simplicity  the  single  matrix  elements  of 
the  form:  a)  •£  b)  |e“*k r|t^>  and 

c)  e”**^|i|^>  separately: 

a)  By  using  commutation  relations  with  the  hamiltonian  of  the 
system*) 


H  -  f  +  V(?) 

1  ’  2_p, 

we  have  ,  ,  ,  .  . 

<Tfl??e’a> Ifi)  =  L  <if'CS  pe 


(4) 


«jm>- 


U-p 


[W 


it? 


i  f{>  mmcs*, 

HJ,  H]l*> 


nm  * 


Now  it  is  clear  how  to  proceed.  To  understand  what  we  have  done. 


*)We  assume  in  this  paper  that  the  potential  energy  V(r) 
does  not  depend  on  velocity  (P/p.)  and  spin  (if)  .  It  is,  however, 
straightforward  to  calculate  the  modifications  of  the  following 
equations  if  this  assumption  is  not  made. 


and 


-ifr 


HJ ,  ft]  *  i  VJe  * 


t«  (&.&)[*■*,  H] 

“  ^  (7) 


^  Jr  /  / 

where  k  =  1^-  is- a  unit  vector.  Equation  (5)  now  r4ads 

p'-^l *  >  =  4  <% IfrP.VJ e*4*/  * >  * 


.  e-'*>/K>  t-  a  /rr  ^  *  &  )^r ' #  > 


(8) 


We  see  that  we  have  an  expansion  of  our  matrix  element 

in  a  power  series  of  7  and  .  The  last  term  in  (8),  in  which  -■ 

c  pc2 

no  potential  occurs,  is  so  to  say  a  "rest  term".  It  is  straight¬ 
forward  to  calculate  all  higher  order  corrections.  The  first  term 
in  (8)  ia  just  the  result  of  Bloch  and  Nordsieck4*5) We  shall 
come  back  to  this  statement  later.  ; 


b)  We  can  .proce.ed  similarly  to  the  outlined  method  in  evaluating 
all  other  types  of  matrix  elements.  We  get 

a/e-^/VJ)  «  l  <%l - 
=  < %ICcr^,  HJ, fill Vs>  -£< % 


Hjpre  we  can  check  once  more  the  earlier  statement  that  the  second 
term  of  (2)  is  at  least  of  the  order  ~  compared  to  the  first  one. 


e) 


!#  it  e'iTrtr;>  -  &  <%*&*■ ’ 


*  £L  <%! e^l ** > 


From  the  last  formula  we  see,  that  if  in  the  nuclear  potential,  V(?) 
no  spin-dependent  part  is  present,  the  commutator  [<T*  R,h]  will 
vanish  and  the  first;  nonvanishing  term  in  the  matrix  element  will 

s$Jl  .--lx 

be  of  the  order  ^  -  hi.  compared  to  the  first  term  of  (8) 

!  O  JJk 

and  therefore  very  small.  All  these  expansions  (8),  (9)  arid  (10) 
of  radiation-matrix  elements  in  power  series  of  the  quantities 
-  and  have  the  advantage  that  the  . nuclear  potential  V(i*) 

C  pt 

occurs  in  the  expanded  matrix  elements.  Therefore  one  can  pro¬ 
ceed  for  short  range  potentials  as  in  usual  nuclear  radiation 
theory  and  a  multipole  expansion  becomes  meaningful:  Only  those 
regions  of  space  where  the  potential  differs  appreciably  from 
zero  will  contribute  to  the  radiation. 

1.2.  Nuclear  Bremsstrahluno-Matrix  Elements  for  Spin  Independent 
Nuclear  Forces. 

In  the  following  we  shall  discuss  a  simple  and  special  model 
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for  nuclear  bremsstrahluang.  It  is  our  aim  to  show  in  this  section 
how  the  formalism  works;  later  in  part  II  of  the  paper  we  shall 
compare  the  effect  of  different  nuclear  potentials  on  the  brerns- 
strahlung.  Here  we  treaf  the  two  particles  without  spin  and 
assume  that  the  internuc Aeon  force  is  static  (not  spin-or  velocity- 
dependent).  As  we  have  ^ust  remarked  following  equation  ( 10 ) a 
the  radiation  according  -fcr o  the  Pauli  term  is  indeed  very  small 
for  such  potentials.  So  we  are  consistent  in  neglecting  the 
Pauli  interaction  with  thx.  e  rad  iation  field.  We  further  like  to 
study  the  pure  nuclear  brr- emsstrahlung ;  i.e.,  we  are  not  interested 
in  the  bremsstrahlung  du&  to  the  coulomb  potential.  We  therefore 
consider  thp  scattering  of  neutral  particles  by  charged  nuclei. 

Of  courses-  the  target  nuc  leus  should  not  be  too  larges  so  that 
the  reboil  will  be  apprec  iable  .  An  examples-for  this  case  would 
be  the  scattering  of  neutrrons  by  protons  or  o-particles.  Further¬ 
more  we  neglect  many-body  effects,  which  may  probably  be  of  interest 
if  the  target  nucleus  is  xrore  complex  and  not  only  a  proton.  The 
■formalism  is  expected  to  Jcecome  more  complex  when  both  particles 
are  chargeds  naturallys  b«ecause  of  interference  between  coulomb 
and  nuclear  bremsstrahlun«g.  The  interaction  (3)  transformed  ac¬ 
cording  to  (8)  and  (9)  now  reads 


Hint  - 


-et  i  j.  r 

C  AH  *w  L  /*> 


(11) 


10 


We  note  that  one  has  to  use  the  exponential  factors  of  equation  (3) 
in  formulae  (8)  and  (9).  So  where  we  have  a  vector  1c  in  (8)  and 

Wl  n,  t-  1  m  • 

(9)  it  should  essentially  be  replaced  by  -jjp  k  *=  k  . 

For  simplicity  we  choose  here  the  internucleon  force  to  be 
a  square-well-interaction.  Later  in  part  11  we  will  discuss 
another  type  of  interaction.  So  we  have 


V  =  ]/{T)  =  [*>■—*  0<= 

V  1  '  L  0  for  ri 


fir  r  fi  flu 

* 


and 


(fra/lV  *  -V*  fx"  ^  ^r*  f  I* 


(12) 


(13) 


where  j  p  are  unit  vectors  in  spherical  'basis6);  Expanding  the 

plane  waves  in  (11)  in  electric  and  magnetic  multipoles,  we  can 

—2- 

restrict  ourself  to  dipoles,  since  ka,*»  10  when  the  energy 

of  the  incoming  particles  is  pf  the  order  3  MeV  (and  therefore 

the  highest  photoenergy  is  -  3  MeV).  We  see  that  because  of 

(13)  the  first  term  and  the  last  term  in  (11)  contribute  only 

electric  dipoles.  Since  the  last  term  is  of  the  order'  ^  ~ 

smaller  than  the  first  one,  it  will  be  neglected*.  Let  us  remark, 

that  in  general  the  only  magnetic  interactions  in  (11)  are  coming 

from  the  second  and  third  term  and  are  of  the  order  ^  smaller 

c  i 


than  the  dominant  electric  terms  as  long  as  the  inter nucleon 
potential  is  spherical  symmetric?  i>«.,  V  *  V(r)  .  This, can 
easily  be  proved  by  writing  down  the  multipole  expansion  for  the 
plane  waves. 


*As  can  be  seen  from  (9)  the  contribution  to  the  radiation 
ebming  from  the  interaction  of  the  center-of -mass-motion  and  the 
relative  motion  vanishes  exactly  in  dipole-approximation. 


■  11 

We  shall  neglect  the  magnetic  dipoles  too,  because  In  Our 

case  they  are  of  the  order  ^  £  JL  (fo r  scattering)  and 

m^  v  ,J  1 

it  c  *  m  (for  rt-a  scattering)  smaller  than  the  electric  term. 
However  in  the  case  of  resonance  scattering  (especially  for 
broad  resonances),  for  instance.  One  should  look  at  this  ap¬ 
proximation  with  cautiipjfr:  if,  for  example,  we  have  a  broad 
p-wave  resonance,  the-  magnetic  contributions  might  be  appreciable* ) . 
Let  us  use  circular  polarized  light-quanta 


=  -  F  up  m  P  ('frifatO)  [r 


(14) 


p-  * ; *  1 


where  Tfp  is  the  circular  unit  vector^)  in  the  system  with 
z-axes  along  the  photon  propagation  vector  k  •  are 

the  angles  specifying  the  k-direction  in  the  laboratory  system. 

Now  we  have  finally 

<Hint)Dipole  ='fp  ^  ^  <15> 


The  Wave  Functions  for  the  Initial  and  Final  States. 

The  Initial  wave  function  characterized  the  incoming  neutron! 
flying  in  the  z-direction  (axis  of  quantization)* 

I 

f.  -  eih*  i-  ('U’i(n)  -ut(«i))  Yte  (*« ,  ftti)- %o  (?<  f) 


yith  the  asymtotic  form 

«*»  ^  fJMl  tlhr 

r 


(16) 


% 


*Furthermore  one  can  see  from  (4)  that  for  static  potentials,  we 
are  discussing,  these  terms  are  of  the  order  "lCa*  smaller  than  the 

dominant  terms i  because  in  the  approximation  we  have  iic.tr t 

'  [H,  e  J  =  0. 
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The  partial  waves  are 

Ut  («)  ~  ^  ‘  *  <{*  V  , 

.«  r.** 

Ve  rtr.0  •  *  e  [  £  fo,  ^  fr/rj  A  r<«C 


where  f#(K 


.)  5  (cc&fti)  jtfcz)  -  (Ki°^) 

i  c°  = 

N  is  the  number  of  partial  waves  with  phase  shifts  if  (K)  differing 
appreciable  from  zero. 

The  final  state  characterizes  neutrons  moving  in  a  direction  kp 


-d 


+  f  (vt(Kt) -uc(*>)) Z  Tel 

f  fao 


(19) 


with  the  asymtotic  condition 


% 


‘fa  fjye)  &~ik*r 


(20) 


We  note,  that  the  correct  final  wave  function  must  have  an  in- 

7 

coming  spherical  wave').  Therefore  we  have 

-Hri1  jt(K,r) 

~Cifc  ({.)  f  W°  £ (^)  it  ~  ^  ^ 

«*»•«'*  e  /  t  s 

(21) 
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whetfe  K£  and  f^K^)  are  similarly  defined  as  KJ  and  fg(K^)  . 

The  matrix  element  of  interest  is 

<%  /.*>-«)  Z  D'n  Xt l rw&z> 

•(lOfttf*)  (Kit)  ~ u* (*■(*■) Ka, (*<*))  YeJe,o)£c(t«, I, (/»»>»>)• 

-  icl  (Kfij  UeJW)  Ye-,  (»,  o)  Z  c ((-!,!,( lQ""')  YejK) 

(22) 

The  C-coeff icients  are  Clebsch-Gordan  coefficients  in  the  nota¬ 
tion  of  Rose^).  We  see  from  (21)  and  (22)  that  only  partial 
waves  with  phase  shifts  different  from  zero  contribute,  as 

expected*).  We  like  to  note  that  we  have  dropped  that  term  of 

■ 1  \ 

thfe  matrix  element  which  contains  plane  waves  on  both  sides, 
because  it  does  not  contribute  to  the  cross  section.  Formula  (22) 
contains  the  correlations  between  the  scattered  particles  and 

x 

the  emitted  photon.  Averaging  over  both  directions  gives 

rtf.i  = 

-  f  («'*,,!,€  I  "Of  YeXteo)  l»:m 
•Ltx  t  L 

C(t-l,l,eiK>oY-  !%'(*+*■)  j 

(23  J 

*)This  follows  in  this  case  from  the  fact  that  ZJt. 

tt-*0 
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Here  we  have  used  orthogonality  of  the  Y^m  and  Dmp  and 

unitarity  of  Clebsch-Gordan  coefficients.  Introducing  the 
2  2 

quantities  AQ  and  by  using  (18)  and  (21)  we  get 

0  G 

At  =  [(<*!&((,) A  (**)*■ t(Ke*)T  i?tn 

-  'Uvo  (Ki&))jL-2o**  (Jt(*f)+  Jen  (&))(<*»  A  (*+)* 

•  ^  ^  (Kf)  )(u'°  (*<)  j't+tCti*')  ~ 

-tafc/*'.')  ^(Xts,)) (je(KfX)  jt*i  (*>«■))  + 

(24) 


ana  . 

&.,(*/  tte„ (n,-a) T ~i<^>  (&<**) 

.  (or,  A  («f)  jdt'fi.)  -2^A(*f)  (*■,*))  ( ■«»  £r, 

(*) *<-,(«$) 6* . 

(25) 


w 

~  2vwv 


Inserting  all  factors  from  (18)  and  (21)  we  get  for  equation  (23) 

^  b-P  l.wHi/L,,  *“  T7,  oc^V 

r 


YJy)~Am/  l£r  . 


Here  we  have  used  the  relation  ' tin  ‘  r®*** »  .  The 

bremsstrahlung  croSjS  section  is  defined  as  the  transition 
probability  between  initial  and  final  state  divided  by  tfye 


(26) 


15 


velocity  of  the  incoming  particle. 

d<r=  X 


(27) 


where 


r-  sS?  |  «&/«-  <r*> 


(«■  g?  1 


(28) 


The  factor  2  in  the  denaity  of  states  p(E£)  is  coming  from 

the  spin.  L  is  the  length  of  the  normalization  cube.  The 

density  of  states  of  the  photon  is  included  in  the  factor 
r  3  o 

"  the  transition  probability  ).  The  sum  over  the 

8ir7ic 

polarization  £  in  (28)  gives  just  the  factor  2  since  (26) 
is  independent  p  (which  characterizes  the  polarization  --  see 
(14)).  Adding  all  factors  we  finally  get  by  integration  over 
dQ^  dQ^  and  using  (26) 

t-r  -  <&(? 

.  c(t+t,  /,  e  looo)1  A?  +  fre-0  Cfa  '<  w)L&] 

(29) 


C  V 

Let  us  now  consider  the  cross  section  (29)  in  the  limit  j  ^  0  $ 

i.e.,  in  the  limit  of  low  energy  scattering  from  a  weak  potential. 

In  this  case  only  s-scattering  occurs  and  the  s-phase  shift  is 
given  by 
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z*l  ^ 
/gr  $««****/  ^  3 


e%  >wv 


(2)0) 


2  E 
where  x  =  = 


Inserting  this  into  the  expressions 
2 


» v2  =  f- 

to  0  fco 

(24)  and  (25)  for  and  and  developing  jQ(x)  88  1 

jj_(x)  ~  ^  ,  nQ(x)  ~  ~  “  f  we  obtain 

Axt  B,x  **  (*tx  *■*•'*'  *-)  ,  *' 1  t; 


Ol) 


So  we  have  in  lowest  order  for  the  square  of  the  matrix  element 
a  constant  and  for  the  power  radiated  via  bremsstrahlung  according 
to  (29)  a  dependence  like  ft  on  energy.  This  is  shown  schemati¬ 
cally  in  figure  1. 


Fig.  1.  Power  spectrum  for  nuclear  bremsstrahlung 
at  very  low  energies  from  weak  potentials. 

1.3.  Approximate  Bremsstrahlung  Cross  Section  for  n-a  Scattering.- 
We  shall  now  present  an  approximative  calculation  for  the 
bremsstrahlung  cross  section  for  neutrons  scattered  by  He4.  The 
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approximation  made  is  that  the  neutron  spin  is  neglected  and  that 
the  wave  functions  for  this  scattering  process  are  assumed  to  be 
those  of  a  square  well  with  parameters  a  =  3.21  ferrni, 

VQ  =  19.65  MeV  according  to  the  work  of  Sack,  Biedenharn  and 
Breit^)o  We  used  further  the  experimental  phase  shifts  of  Clementel 
and  Villi'1'0),  where  we  identified  their  P^^-inverted-doublet,  phase 
shift  with  our  p-phase  shift.  The  energy  E.  of  the  incoming  neu¬ 
trons  is  chosen  to  be  E^  a  1.4  MeV,  just  above  the  P3/2  resonance. 
The  numerical  result  for  the  energy-dependent  terms  of  the  power 
radiated  by  bremsstrahlung  is  shown  in  figure  2. 

_2 

Using  formula  (29)  one  gets  a  cross  section  of  the  order  10 
millibarn  if  the  energy  resolution  is  chosen  dE£  -  0.2  MeV.  It 
is  interesting  that  the  bremsstrahlung  differential  cross  section 
shows  a  maximum  in  the  region  of  the  p-resonance.  This  can  be  ex¬ 
pected.  However,  the  details  should  not  be  taken  too  seriously, 
because  the  approximations  made  are  crude. 


II. 

Comparison  of  the  Effects  of  Hard-Core  and  Non-Hard-Core  Potentials 
on  Nuclear  Bremsstrahlung. 

We  now  discuss  the  possibility  of  differences  in  the  brems¬ 
strahlung  spectrum  due  to  hard-core  and  non-hard-core  potentials. 
The  formalism  developed  in  section  I  where  the  nucleon  potential 
is  explicitly  worked  into  the  radiation  matrix  elements  suggests 
such  a  possibility.  It  turns  out  that  bremsstrahlung  cross  sec¬ 
tion  is  different  by  only  a  few  [5%)  percent  for  both  types  of 
potentials. 
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II. 1.  The  Wave  Function  of  a  Hard-Core  Potential. 

The  type  of  hard-core  potential  we  shall  investigate  is  $hown 
in  Fig.  3.  The  hard-core  is  represented  by  the  potential^  . 


Fig.  3.  The  type  of  hard-core  potential  which  we 
investigate.  We  take,  however,  the  limit 
Let  us  construct  the  wave  function  of  this  nuclear  potential  in 
the  usual  way  by  constructing  the  wave  functions  in  regions  I, 
II,  III  (fig.  3)  and  fitting  their  magnitudes  and  derivatives  at 
the  boundaries.  We  introduce  the  following  notations 


,  _  i/_v(V„-eT 

I  ^  ■  I  - 


/ 


v= 


(32) 


For  positive  Vv  and  small  or  negative  E  the  a  v  become  complex. 

In  the  other  cases  they  are  real.  The  equation  for  the  radial 

wave  function  R  is 
e 


d'*€  ,  i.  ,  P-  eJW]Rt  -  0 

f  Ctf  +L  ^  J  (33) 


For  the  regions  shown  in  figure  3  we  have  the  following  solutions 
of  (33)  with  regular  behavior  at  the  origins 
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i«  At  ^WfO 

III  St  (•*,  (h)  -  ^  ^  (34) 

m,  Ct  (<*»  A  '!V*  (ft)  3  h>  **♦  <&>-) 

J  and  n.  are  the  Bessel  and  Neumann  functions  with  the  asym- 
6  6  * 

totic  behavior^) 

-  f_g  — 

+*\'"Uen) 

(35) 

- s . t  •"  ac "0 

p  C+l 

j~  Cf-  i(t+0TrJ 

(36) 

1  *vU  [f>-£(C+')7r] 

The  2/^  are  the  phase  shifts  of  the  e’th  partial  wave  in  the 
regions  indicated  by  the  index  V  •  From  (34)  we  get  the  following 
equations  for  determining  the  coefficient: 
the  continuity  of  the  wave  function 


B  and  due  to 

e  e 

(37) 


i*(  f) 

f  ">0 

ne(f) 

and 

«t*  <f) 

^  ->«0 

and  for  the  logarithmic  derivatives 
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*1  faxA)  *'t(*l*  )1  Z^?  *VW^’« 

•»«%  l _,Xl«*V  JiM*iCcsA)]  £*>i?e"£M  } 


Analogous  equations  hold  for  the  boundary  i/ll.  We  have  there 

^  J.&  fa  6~)  Se.  **•«  ^  J 


and 


^  <f< i (r )  _  ^  £c*o$t.  ^  ^ /V*.  ^  v£  ^efa^ll 

6-)  [ooo  j!e  fax  t)  ~  a***  i fo,(£l/ 

From  the  last  equation  we  can  determine  the  phase  shift 


**»- 


<*J 


Je  fa'lfi')  M/tfaG')  —  J'tfa  &) 


(39) 


(40) 


(41) 


Since  we  are  interested  in  the  limiting  case  V*  -*°° }  i.e., 
infinite  hard  core,  we  have  to  look  at  equation  (41)  in  the 
limit  -*  i°o  .  Using  (36)  we  get 

r  -  fa l  -  K  (CH)Tr )  coo  fa  £-%(t*)r)\ 


itM^} 
fa  tr) 


ioO 


C^'i  C-  t~  %.  (z+i  )  y ) 

»<1  0- 


“  -  (o^o  l  -  %  i^Orr)  -  ^ 
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Remembering  that  a-^  is  complex  for  positive  ,  we  put 


(42) 


Inserting  this  in  the  arguments  of  and  writing  this  in 

exponential  factors  we  find 


and  therefore  we  have 


<•**>  £  **{<*) -tts fit*) 


j-cfai.t')  Ei  __  *cx  toe'kif/Sl 

KjiUii)  L  '  JtM)  HcfejV'-I 


(44) 


During  the  last  step  we  kept  the  range  of  the  hard  core,  "b*1, 
fixed.  From  equation  (38)  we  can  calculate  the  other  phase 
shift  ^  .  We  have 

■  n(j)_  jfc  UjA )  -  Te  |e.  4 ) 

^  &  Hi  -  re  n& 

*3  0«  ~  ^  («zA)] 


(45) 


(46 ) 


where 
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11*2.  The  Wave  Function  of  the  Square-Well  Potential. 

These  wave  functions  are  well  known.  Although  we  have 


already  used  them  in  section  I8  we  shall  summarize  the  results 


in  order  to  be  complete.  According  to 
figure  4  there  are  two  regions  which  we 
call  IX  and  III.  The  wave  functions  and 
phase  shifts  are  evaluated  with  the  same 
method  as  in  II  1. 


Fig.  4.  Square  well 
potential.  The  parameters 
a,  V*  are  somewhat  dif¬ 
ferent  from  those  of  the 
square  well  in  Fig.  3. 


We  have 

nt  Be.  JtC^r) 


iiu 


c*  ( 


■ 


j-e  ^1°'  )  it-  (<1®,) 

- — - - - 

5l  LcC^i  o£)  ~  </e  (*i&)  toe (*i  3-) 

C  OVr*  (<S&)  ~  ***1  (*3  * 

(<*!.(*,  ) 


(48) 


(49) 


II. 3:  Calculation  of  Bremsstrahlunq  Matrix  Elements  for  the  Two 
Kinds  of  Potentials. 

As  outlined  in  section  I  we  used  the  gradient  of  the  potential 
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for  calculating  the  dominant  terms  of  the  bremsstrahlung  matrix 
elements.  We  have  for  hard  core  and  pure  square  well  respectively. 


i.  'ffcul  1/  =■  F  [  IA,  +-14  i^(V  -a.)  J  Z- 

2:  c^cLV  -  F  ^‘n  (?)  'Kf* 


where 


P=-p^  "t£t  (compare  equation  (15)), 


p-c  ^ 

The  matrix  elements  we  need  are  of  the  type 

<%($)  ly^oU/IVc  (e;)> 


(51) 


where  ^(E^)  and  ^(E^)  are  the  final  and  initial  wave  func¬ 
tions  we  have  constructed  in  equations  (21)  and  (18)  but  with 
some  modifications  in  the  radial  part  of  the  partial  waves 
according  to  the  results  of  II. 1  and  II.2.  We  further  make 
some  changes  in  the  calculation  of  the  bremsstrahlung  cross 
section:  We  neglect  the  contributions  of  higher  partial  waves 
(which  we  included  exactly  in  section  I),  Therefore  we  write 
for  the  wave  functions 

f'(e)  =  £  ve.  (e) 


This  simplifies  the  calculations  in  the  hard  core  case*  which  we 
now  only  do  and  we  note,  that  the  square  well  is  discussed  in 
section  I.  Since  the  angular  dependent  terms  are  the  same  in 
both  cases  (50)  we  only  calculate  the  radial  parts  of  the  matrix 
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elements.  We  need  the  wave  function  at  the  points  T  =  b,  A, 
as  can  be  seen  from  (50).  At  T  =  b  . we  have  according  to  (34) 

Be.  (V«  ^  ^  ^iXt £))  ^  _> ^ 

B*  C,c*Ve  (  5J  t——}  +  ^  rt  /( 


/  *  ^ 

s  the  Vg. 


We  used  (44)  in  getting  this  result,  vt  i 
So  we  get  for  the  radial  part  of  the  matrix  element  between 
the  partial  waves  (E)  ,  %  , (E.  ) 


c  c  •*•  /  \2_  ♦  '2 ?[(» 

a IE)  1 14  *'-6)  /  (f. ;^w=  4 

y_7»Wi  ,  A  *!*££* 

(54) 

Inserting  the  factors  a2  ,  from  (32)  and  (42)  and  the 

coefficients  B*  ,  B  ,  from  (47)  we  get 

<%(e)lK.f(r-t)IWn)>~i,f*4-  FKtFfWi  ct  4' * 
x  Qt  («,A)-tj  He  (<}*))  (je  W*  )  -  i  A  'j«<  • 

.  .  (55) 

*  Ut(£)  Uc‘(e0  '  o 


where 


/-_  j*'M)  7 

.  //r)  _  *z  2e 


(56) 
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All  quantities  in  (55)  depend  on  energies  E  or  E.  .  The  C 
can  be  taken  from  (21).  The  other  contribution  to  the  matrix 
element  is  M  .  "N 


CMV  I  K  %'  =  V%  C*Cl‘ 

.Qti  (<,  A  -  ^  j h'(i>  ht'faiA))  to*  'A 


(57) 


Here  we  have  used  the  last  of  equations  (34).  By  adding  (55) 
and  (57)  we  get  finally  for  the  radial  part  (of  the  hard  core 
plus  square  well  potential  contribution  to  the  matrix  element). 


(WM-h  ^Vmj) 


We  see  here,  that  the  only  essential  difference  of  thebrems- 
strahlung  matrix  element  of  case  (50.1)  (hard  core  square  well) 
compared  with  case  (50.2)  (pure  square  well)  is  contained  in 
the  first  term  of  the  bracket  [  ]  in  (58).  The  functions 
(56)  describe  the  contribution  to  bremsstrahlung  from  the  hard 
core.  It  is  interesting  to  note  that  the  contributions  of  the 
hard  core  and  the  square  well  in  (58)  are  opposite  in  phase; 
i.e.,  they  interfere  destructively.  Of  course,  the  square  well 
is  accelerating  the  particles  in  the  opposite  direction  to  that 
of  the  hard  core. 

It  is  now  easy  to  write  down  the  radial  part  of  the  same 
matrix  element  in  the  case  (50.2).  This  we  do  for  completeness: 
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<r  %  m  i  K  ter-*)  i  %>  (?<■)  -  <r<v  fr*  &««-***•) 

«(y  65fi ) -  Vj  T^e-Vt*  c«  /e w«^>  (!*> 


II. 4.  Fitting  the  Parameters  of  the  Potential. 

The  range  of  the  hard-core  potential  "b"  is  taken  from 
high  energy  scattering  data  )  to  be  b  -  0.4  fermi.  The 
parameters  a  ,  and  a  ,  are  in  both  cases  somewhat 
different.  Let  us  assume  a  ,  to  be  given  (for  instance 
from  low  energy  scattering  data  by  fitting  the  scattering  length 
and  effective  range),  we  then  determine  the  parameters  a  , 
so  that  we  get  for  the  hard  core  the  same  scattering  length  and 
effective  range  as  for  the  square  well  potential.  In  other  words, 
we  expand  the  s-phase  shift  in  a  power  series  of  the  energy  and 
require  that  the  two  lowest  terms  of  the  series  are  the  same 
for  both  potentials.  We  then  get  two  equations  which  determine 
the  parameters  a  ,  in  terms  of  the  parameters  a  ,  V2  . 

The  phase  shifts  are  given  by  equations  (45)  and  (48);  however, 
the  s-phase  shifts  in  both  cases  are  simply  connected  by 


Let  us  now  introduce  some  new  quantities 

=  Vvi+*1’  = 

=  X  1  °<ib  s  x,  *  ff;  i  (61> 
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The  analog  definitions  hold  for  the  barred  quantities.  So 
we  have  expressed  the  old  quantities  a  ,  V£  ,  E  in  terms  of 
the  quantities  SQ  ,  y  ,  x  .  The  equation  for  determining  the 
unbarred  quantities  from  the  barred  ones  reads 

^ Cu ,  y/x)  =■  ofy  y/x)  =  ^  9) 


(62) 


Let  us  now  be  more  specific  and  restrict  ourselves  to  the 
case  of  neutron-proton  scattering.  Here  the  parameters  are  due 
to  the  results  of  Blatt  and  Jackson"^)  (1950):  a  =  2.58  f  and 
^2  =  -  13.3  MeV  .  Therefore  we  get  **  5.67  MeV,  y  w  1.531  , 
ss  41.602  .  Keeping  these  quantities  in  mind  we  proceed 

and  get  from  (62) 

(i) 


We  now  can  easily  derive  from  (48)  that 


(63) 


fi 


(64) 


Since  we  expect  y  s  y  ,  dfagy  *  ctgy  =  0.04  and  in  the  low  energy 
limit  we  have  ^  «  1  .  So  we  can  develop  expression  (64)  and  get 


(65) 


Now 


we  develop  (63)  also  'since  «  1  and  get 


~  0*4 $("&)>&]  (66) 
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We  therefore  can  write  equation  ((?2)  in  the  form 


This  equation  determines  the  parameters  y  , from  the 
parameters  y  ,  .  We  now  assume  that  the  square  well  does 

not  differ  much  in  both  cases;  i.e.,  we  assume  that 


14=  14.  (4+  f) 
So  -f.  ('I **}) 


46.  'I 


(68) 


are  to  be  determined  and  are  assumed  to  be  small  numbers. 


^  and  J 

It  turns  out  that  these  assumptions  are  consistent  with  the  result 


we  will  get  for  ^  and  I j  .  We  now  remember  the  definitions  (61) 
and  introduce  (68)  in  (67).  Setting  the  first  two  coefficients 
of  the  same  power  in  energy  E  in  (67)  equal  and  linearizing 
the  two  equations  for  and  yj  we  get 


^ii  f  +•  dixy  * 

where  (69) 

a„  -  =  I*  ' 

c„~o  ,  c*  =  /?* 
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From  these  equations  it  is  easy  to  determine 
we  get  |  =  0.129  ,  fj  =  0.125  which  is  consistent  with  our 
assumption  that  these  quantities  should  be  small.  According 
to  (68)  we  have  the  result  that  the  two  square  wells  are  nearly 
equal*  The  square  well  connected  with  the  hard  core  is  somewhat 
deeper  and  somewhat  less  extended  than  the  other  one.  Further 

o  —  — o 

it  is  easy  to  see  that  V2a  =  ^2a  to  a  verV  900C*  approximation. 

So  we  get  the  same  binding  (or  virtual  binding)  energy  in  both  cases 
What  we  have  in  fact  done  is  to  demand  equal  scattering  length.  a,nd 
effective  range  for  the  two  potentials  we  are  discussing. 


f 


and 


f 


arid 
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RESULTS  AND  DISCUSSION 

We  have  done  some  numerical  work  in  order  to  find  out  differences 
in  the  cross  section  for  bremsstrahlung  due  to  the  two  types  of  po¬ 
tential  we  are  discussing.  The  result  is  shown  in  Figure  5.  The 
two  potentials  give  for  s-bremsstrahlung  nearly  the  same  shape,  and 
the  cross  sections  seem  to  be  different  in  both  cases  by  only  a  few 
percent  The  bremsstrahlung  of  the  hard  core  potential  is 

smaller  than  that  of  the  pure  square  well.  This  is  due  to  the  de¬ 
structive  interference  of  both  potentials.  The  particles  are  accele¬ 
rated  by  the  hard  core  in  the  opposite  direction  as  they  were  accele¬ 
rated  by  the  square  well.  As  can  be  seen  from  (58)  the  shape  of  the 
bremsstrahlung  spectrum  is  for  both  cases  nearly  the  same,  because 
the  shape  is  essentially  determined  by  the  squared  factor  before  the 
bracket  [  ]  and  this  is  in  both  cases  nearly  the  same.  The  bracket 
l  ]  depends  only  smoothly  on  energy  and  is  nearly  constant  over  the 
energy  interval  we  have  considered.  However,  it  seems  that  the 
nuclear  bremsstrahlung  should  be  measurable,  since,  with  an  energy 

resolution  of  about  dE^  =  0.2  MeV,  we  have  for  the  differential 

-2  -3 

cross  section  do  *  10  -  10  millibarns*. 

At  low  energies,  before  the  cross  section  drops  to  zero,  it 
shows  a  maximum  (E£  «  50  KeV).  This  is  due  to  the  bigness  of 
the  s-wave  function  at  the  nucleus  (virtual  bound  state)  and 


*It  seems  useful  to  note  here  that  the  magnitude  of  the  cross 
section  depends  on  the  energy  of  the  incoming  particle.  So  per¬ 
forming  experiments  and  comparing  these  with  theory,  one  has  to 
know  several  experimental  parameters  accurately. 
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similar  to  the— maximum  of  the  magnetic  photo-disintegration  cross 
section  at  the  same  magnitude  of  energy. 

We  have  done  the  calculations  without  spin.  This  gives  good 
results  as  long  as  we  discuss  static  potentials;  i«e.,  without 
spin-  and  velocity-dependent  terms.  This  statement  has  been 
proved  following  equation  (10).  However  in  the  next  step  one 
should  discuss  this  other  class  of  nuclear  potentials  and  try 
to  find  out  possible  differences  in  the  nuclear  bremsstrahlung. 
Further  it  seems  to  be  an  advantage  to  look  on  bremsstrahlung  at 
even  higher  energies,  because  one  can  expect  there  a  bigger  cross 
section  since  more  partial  waves  contribute.  This,  however,  is 
more  or  less  only  a  numerical  task. 
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